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In light matter interfaces based on the Faraday effect quite a number of quantum information 
protocols have been successfully demonstrated. In order to further increase the performance and 
fidelities achieved in these protocols a deeper understanding of the relevant noise and decoherence 
processes needs to be gained. In this article we provide for the first time a complete description of 
the decoherence from spontaneous emission. We derive from first principles the effects of photons 
being spontaneously emitted into unobserved modes. Our results relate the resulting decay and noise 
terms in effective equations of motion for collective atomic spins and the forward propagating light 
modes to the full atomic level structure. We illustrate and apply our results to the case of a quantum 
memory protocol. Our results can be applied to any Alkali atoms, and the general approach taken 
in this article can be applied to light matter interfaces and quantum memories based on different 
mechanisms. 
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I. INTRODUCTION 

The strong and coherent interaction of light with mat- 
ter is a prerequisite for many approaches towards quan- 
tum information technologies. In particular long distance 
quantum communication relies on efficient light-matter 
interfaces which allow for a coherent transfer of quantum 
information from light to stationary carriers and back 
[1]. Also architectures for quantum computations based 
on light will depend on efficient light-matter interfaces 
for buffering and storing quantum information carried 
by light. 

Optically dense atomic ensembles have proven to be a 
particularly promising technique for achieving strong co- 
herent light matter interaction. Thereby quantum states 
of propagating pulses of light are mapped onto states 
of collective atomic (pseudo) spins. This essentially re- 
quires a mechanism to fiducially and reversibly convert 
photons into ground state spin excitations with long co- 
herence times. Several mechanisms have been explored 
for realizing such an atomic ensemble based light matter 
interface, e.g. Raman transitions, Electromagnetically 
Induced Transparency [2-4], Spin Echoes [5], and Fara- 
day rotation, see [6-8] for comprehensive reviews. The 
Faraday effect — which will be the topic of this article — 
consists in the rotation of light polarization depending on 
atomic spin polarization, and vice versa. In the context 
of light matter interface it has been successfully used to 
create squeezed states for spin [9, 10] and light [11, 12], 
entangled states of collective atomic ensembles [13, 14], 
quantum teleportation from states of light to atoms [15], 
and quantum memory for light [16, 17]. Apart from the 
light-matter interface the Faraday interaction has impor- 
tant applications also in continuous nondemolition mea- 
surement of atomic spin ensembles [18, 19], quantum- 



state control/tomography [20], and magnetometry [21- 
23]. 

These protocols can all be understood in terms of a 
rather simple model of the Faraday interaction. In this 
model atoms are assumed to have a spin 1/2 ground state 
and a spin 1/2 excited state [24]. Far off resonant light 
probing this dipole allowed transition will then experi- 
ence a polarization rotation due to the dipole selection 
rules for such a 1/2 — 1/2 transition. By the same effect 
the atomic polarization will be rotated by light. For far 
off resonant light atoms will be only very weakly excited 
such that this birefringence of the atomic medium can 
be understood as being due to the polarizability of the 
atomic ground states alone [25, 26]. The coherent inter- 
action of light with atoms arises from the real part of the 
atomic polarizability. By the Kramers-Kronig relation it 
is clear that the corresponding imaginary part will neces- 
sarily be non-zero and add some incoherent effects to the 
dynamics. Sure enough, these effects can be understood 
as resulting from spontaneous emission events. They will 
cause both, decoherence of light (absorption) and deco- 
herence of atoms (spin decay) . While both effects can be 
kept small as compared to the coherent dynamics, they 
are ultimately unavoidable on a fundamental level. 

Both decoherence effects — light losses and spin de- 
cay — can be included on the level of the simple spin 1/2 
model giving a qualitative understanding of the trade- 
off between coherent and incoherent contributions to the 
light matter interaction [27-29] . However, in view of the 
experimental achievements it will become increasingly 
important to further gain a detailed and more quanti- 
tative understanding of these effects and the resulting 
tradeoff. A theoretical description based on a realistic 
atomic level scheme, such as the one shown in Fig. 1, 
and starting from first principles has not been given so 
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FIG. 1. Level structure of Cesium. We assume that our laser 
light is blue detuned by an amount —A from the transition 

f = 4 ^ F' = 5 



far. In this article we will provide such a description. 

From the standard dipole interaction of a single multi- 
level atom with the three dimensional electromagnetic 
field we consistently derive effective equations of mo- 
tion for the collective ground state spin and the forward- 
propagating light modes relevant to the description of 
the light matter interface. In the derivation we will keep 
track of the effects of events where photons are emitted 
to any other than the forward direction. This eventually 
adds decay and noise terms to the equations of motions 
whose origin and dependence on the details of the atomic 
level structure are fully understood and explained for the 
first time in the present article. Important steps towards 
a deeper understanding have been taken before in [30]. 
The insight gained by the detailed knowledge of losses 
and decoherence will enable an optimized operation of 
the light matter interface based on the Faraday interac- 
tion, and therefore contribute to boost its performance 
and fidelity. 

The article is organized as follows: In section II we 
will briefly summarize the Faraday rotation effect and 
connect it to the atomic polarizability. In section III we 
give a description of the coherent part of the dynam- 
ics, deriving effective equations of motion for collective 
atomic spins and forward propagating modes. We will 
write down and solve these equations of motion also for 
canonical operators for light and atoms, and illustrate 
the working principle of a quantum memory protocol. In 
section IV we arrive at the main results of this article 
and include in our derivation spontaneous emission and 
decoherence starting from first principles. The result- 
ing equations of motion can again be expressed and inte- 
grated in terms of canonical operators. As an illustration 
we apply the result to the quantum memory protocol. In 
section V we provide a self-contained executive summary 
of our results. Readers just interested in applying the 



correct model for decoherence in the Faraday based light 
matter interface can directly consult this section. 



II. EFFECTIVE INTERACTION 

We consider alkali atoms, which have a single electron 
outside a closed shell. This electron will be in a 51/2- 
state, giving two stable ground levels: a higher - and a 
lower hyperfine manifold, describing respectively a state 
where the electronic spin J = ^ (in this article we put 
h ^ c ^ 1) is parallel to the nuclear spin / giving a total 
oi F = I + ^ and a state where the spins are pointing 



opposite giving F = I 



see Fig. 1 where we show 



the level structure of the specific case of Cs where / — 
I so that F = 3, 4. There will be two dipole- allowed 
transitions to the excited states P1/2 and P3/2, each of 
which consists of several hyperfine states with spin F' 
(primed variables always refer to electronically excited 
states). In this article we will mainly consider the case 
where atoms are initially prepared in the Si/2{F = / + 
|)-state and are driven on the D2-transition to P3/2 as 
indicated in Fig. 1 for the example of ^'^^Cs. 

In the standard dipole and rotating wave approxima- 
tion the coupling of light to atoms is given by 



_(d(+)E(-)+d(-)EW^ 



where E^^^ and d^*' are the positive/negative frequency 
component of, respectively, the electric field E = E+ -|- 
E^ and the dipole moment operator d = d+ -I- d~ . We 
will be interested here in the interaction with far off- 
resonant light only. In this case the excited states will be 
only very weakly populated. In the limit of low satura- 
tion excited states can be adiabatically eliminated (Ap- 
pendix A), and the light-atom interaction is described by 
the effective Hamiltonian 



(1) 



We have here introduced the polarizability tensor opera- 
tor 



a. = 



E 



PgdPp'dPg 



where the projection operators are defined as 

PF=Y,\F,m){F,ml Pg = Y,PF, 

rn F 

Pp. = ^ m') {F\ m'l , Fe - E 



(2) 

(3) 
(4) 



such that Pg + Pc = 1. Here Ap' is the detuning of the 
light from the S1/2F — )• P3/2F' transition, see Fig. 1. 
The polarizability operator a in (2) has to be under- 
stood as a 3 X 3 matrix whose Cartesian components are 
given by a.ij = — X^f' Pg^iPp'^jPg / Ap' ■ Each of these 
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components is an operator acting exclusively within the 
subspace of ground states Si/2{F = I ± 

The effective interaction in (1) describes second or- 
der processes where a photon is absorbed and reemitted, 
while the atom makes a transition from a ground to an 
excited state and back again to a ground state. In such 
a process the initial and final ground states can cither 
have the same or different total spin F = / ± ^. In 
order to describe processes involving transitions within 
the F-subspace we introduce Hpp = PpHf^^Pp. Tran- 
sitions between F o F — 1 we describe by HpF-i = 
Pp-iH^^^Pp + PpHf^^Pp-i- In general such second or- 
der processes could also give rise to changes of F by 2, 
but since we are considering 6*1/2 ground levels there are 
only two ground levels, e.g. as shown in Fig. 1 for Cs 
F = 3,4. 

Consider first transitions within one F-subspace as de- 
scribed by the Hamiltonian HpF- Due to the fact that 
we have conservation of angular momentum in the in- 
teraction, we can decompose the polarizability operator 
apF = PfolPf into its tensor parts [20, 26], namely 



apF = - ^(ao + iai'i x -|-a2Q). (5) 

Here the dipole matrix element is defined as dp = (2 J' -|- 
l)|(J'||d|| J)p with J and J' being the electronic angular 
momenta of the ground and exited states, respectively 
(see Appendix B), A is the detuning from resonance as 
shown in Fig. 1. The three terms are the scalar, vec- 
tor and second rank tensor part, respectively, j is the 
spin operator for the total ground state spin, that is 
j^|F, m) — F{F + 1)|F, m). In the vector component 
we use here the short hand notation 



jx 



which is equivalent to taking the cross product of j with 
the vector to the right and then the scalar product with 
the vector to the left. One can use vector product prop- 
erties in order to transform the expression as follows 
iE(-) • [j X E(+)] = -ii ■ [E(-) X E(+)]. The latter form 
could be more convenient since it is just a scalar prod- 
uct of the spin and a vector for light which resembles the 
Stokes operator characterizing the polarization of light, 
as will be discussed in more detail below. The second 
rank tensor part Q is defined component-wise as 

2 

Qij = -Ud] +3]ii) + 

The tensor decomposition of ol is given in detail in Ap- 
pendix B where we also give the explicit expressions for 
the real coefficients ag , ai and 02 and their dependence 
on the detuning. In Fig. 2 we show these coefficients for 
the case of Cs. Note in particular that the second rank 
tensor polarizability vanishes for large detunings. 
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FIG. 2. a and b coefficients as functions of detuning —A 
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For the other case where the final spin state has a 
different total spin than the initial state, we have 



a™ = -^(6iTW+62T(2)) 



(6) 



The coefficients 61 and 62 and the irreducible tensor op- 
erators T^^' and T^^^ are given in Appendix B. We do 
not have a scalar tensor component here because in these 
processes the spin state is changed such that the corre- 
sponding Hamiltonian cannot have a component propor- 
tional to the identity operator i.e. a scalar component. 
In Fig. 2 we show bi and 62 a function of the detuning. 
Note that hi and &2 have the same scale as ai and 02, 
respectively. 

The effective interaction Hamiltonian (1) describes 
elastic Rayleigh and inelastic Raman scattering of light 
on atoms. It is a Hermitian operator and as such gives 
rise to a fully coherent evolution for the overall system 
comprised of the atom and the quantized electromagnetic 
field. The quantum coherent effects resulting from this 
effective interaction have been widely used in the field of 
quantum information, both theoretical and experimental 
[16]. Note however, that within the regime of its valid- 
ity — that is far below saturation of excited states — 
the interaction Hamiltonian (1) still contains and cor- 
rectly describes the effects of spontaneous emission on 
both atom and light: If light is treated as a reservoir and 
we trace over the random position of the atoms we will 
show that we can describe spontaneous emission which 
leads to decoherence effects within the ground state lev- 
els. Vice versa, the atom (or an ensemble of atoms) can 
provide an effective absorptive medium for light. In the 
following sections we will first treat the coherent inter- 
action and then include the decoherence effects resulting 
from the effective interaction Hamiltonian (1). 
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FIG. 3. The spin polarized atomic sample interacts with the 
off-resonant light field Ep and non-forward modes Es in vac- 
uum states. 



III. COHERENT INTERACTION 

A. Interaction with a single atom 

We consider now a situation as shown in Fig. 3. A 
cloud of atoms interacts with an incoming pulse of light 
described by the forward propagating electric field Ep, 
and we are interested in the coherent evolution of the 
spin of atoms and the forward propagating field modes 
[31, 32]. The total electric field E = Ep Es will be the 
sum of forward propagating modes and the remaining, 
non- forward propagating modes Eg, cf. Fig. 3. These 
modes will be treated as a reservoir in the later sections 
and give rise to the decoherence effects to be discussed 
there. 

For the moment we will thus restrict ourselves to light 
propagating along the z-direction. In principle the for- 
ward modes will also experience some spreading [33, 34], 
but here we take the opening angle to be roughly zero, 
and use a one dimensional description of Ep(z). The sep- 
aration into a forward propagating mode Ep described by 
a one dimensional equation of motion and a set of non 
forward modes Es is further justified in Ref. [6, 30] and is 
valid if the Fresnel number of the ensemble is much larger 
than unity. In this model the field is described in terms 
of position dependent annihilation (creation) operators, 
defined by 

a.(z,t)= J ^afe.(^)e^'=^ (7) 

where a = x,y labels the transverse polarizations. These 
operators obey [aa{z,t),al.,{z' ,t)] = Scrcr'S{z — z'). In 
terms of a^iz, t) and ay{z^ t) the Stokes operators S'i(z, t) 
can be introduced 

Sz^jXalay-ala^), S'o = ^(4a^ -I- ^a^), (9) 

which are a convenient tool to describe the polarization 
state of light propagating along the z-direction. They 
obey the commutation relations for an angular momen- 
tum density [5*4 (z), Sj{z')] = i5{z - z')eijkSk{z). 



Each part in the decomposition (5) of the atomic polar- 
izability operator will give rise to a term in the effective 
interaction Hamiltonian, such that Hpp = i/'"^ +H^^'' + 
^(2) corresponding to the scalar, vector and tensor parts 
respectively. These terms can be conveniently expressed 
and interpreted in terms of Stokes operators. The scalar 
Hamiltonian — the tensor-0 term is 

= gaoSo. (10) 

Here we have defined q — — where A is the de- 

tuning from resonance as shown in Fig. 1, uiq the atomic 
transition frequency, and A the beam cross section. It 
can be interpreted as a Stark shift, which equally shifts 
all atomic energy levels proportional to the light inten- 
sity. Conversely it can also be interpreted as an equal 
shift of the frequency for all light modes (off-resonant 
with a detuning A), that is as a new index of refraction 
seen by light. 

The vector Hamiltonian, which for our purpose is the 
most interesting part, is 

= gaiSJ,. (11) 

Sz describes the circularity of light and jz the z- 
component of the atomic spin. In the interaction the 
atomic spin is rotated around the z-axis, by an amount 
proportional to Sz. Likewise, the Stokes vector is rotated 
about z-axis, by an amount proportional to jz- This is a 
circular birefringence effect and this interaction gives us 
the desired Faraday interaction. 

Finally we also have a complicated tensor Hamiltonian 

^ ~ga2iS,i£-fy) + Sy{j,,jy} + 2Soi3fz-f)m. 

(12) 

This amounts to a dynamical Stark shift. The effect van- 
ishes for large detunings since the coefficient 02 goes to 
zero as shown in Fig. 2. It can be interesting for coher- 
ent dynamics in some experiments for tomography and 
collective squeezing purposes [14, 21, 35-37]. We keep 
it for calculation of spontaneous emission, and neglect it 
for simplicity in the coherent dynamics. 



B. Interaction with an atomic ensemble 

So far we have considered a single atom. We will now 
assume that light is interacting with a cloud of Na atoms 
of length L and constant atomic density p as shown in 
Fig. 3. For the atomic spin we therefore introduce the 
continuous spin density operators, which are a sum of the 
single atom (j") angular momentum operators evaluated 
at their respective positions 

jk{z) = '^S{z - Za)jk, k^x,y,z. (13) 

a=l 
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Then the commutator is 

Na 

j„(z')] = ^ ^(Z - Za)S{z' - Zb)[j^, j^] 

= ■iemnkjk{z)S{z - z'), (14) 

in perfect analogy to the commutation relation for the 
Stokes operators. For details regarding the definition of 
these continuous spin operators, in particular concerning 
transverse effects, we refer to [6, 30]. 

We are now in a position to describe the coherent inter- 
action of the forward propagating light modes with the 
atomic ensemble. We ignore here the contributions from 
since it merely accounts for the overall refractive 
index. The effective interaction can then be written 

^fh = 5 r(7-S)(z,t)dz, (15) 
Jo 

where 7 = {-a2{jl - jl),-a2{jx,jy},aijz) is a vector 
describing atomic polarization. In Appendix C we show 
how the time derivative of the light operators can be 
transformed to a derivative in position (z), such that 
when we form the Heisenberg equation of motion we get 
from the relation [Si{z), Sj{z')] ~ i5{z — z')€ijkSk{z) that 

£s(z,t)=g(7xS)(z,t). (16) 

If we replace the operators in 7 with their expectation 
values, we can see that in the interaction the Stokes op- 
erator S gets rotated about the vector 7. Written out in 



full detail the equation reads 

f) (Sx\ ( -aij:, ~a2{jx,3y}\ (Sx\ 

^\Szl \a2{3x,Jy)^a2{]l^jl) j\Szj 

(17) 

These equations have first been derived by B. Julsgaard 
[31]. The rotation of S can be seen to be composed of a 
big rotation (proportional to ai) around the z-axis and 
proportional to the atomic spin along z and a small rota- 
tion (proportional to 02) in the (x, y) plane by an angle 
which depends on the relative angle between the mean 
atomic spin and the Stokes vector. If we only consider the 
oi-terms above, and assume the x component of Stokes 
vector to be much larger than the other two projections 
(as is the case for x polarized light), and assuming the 
rotation angle to be small, then we arrive at the effectice 
equation of motion 

J^g)(M) = .<..&(^j). (18) 

These simpler equations tell us that Sz is conserved in 
the interaction, while Sy receives some contribution from 
the spin component jz- 

A similar analysis can be performed for the the coher- 
ent dynamics of the atoms. We will write the coherent 
Hamiltonian as 

i7coh-g[7-S + 7o5o], (19) 
where 7 is given as above and 70 = —a2j1- From this we 
can determine the coherent evolution of the spin vector j 




a2{3y,iz} -a2{jx,jz}-aijy 
a2{jx,jz} a2{jy,jz} aiix 
-2a2{jx,Jy}2a2{jl-j^) 




-{jyjz}^ 
{jxjz} I Sq 




While the overall form of these equations is rather com- 
plicated we can get valuable insight by treating the differ- 
ent terms separately [31]. For instance if we only consider 
ai-terms as it was done with light and also assume small 
Faraday rotations of atomic spins with strong polariza- 
tion along X then we get the equation 



= (21) 



in perfect analogy to the effective equations of motion for 
light found above, cf. (17) and (18). jz is unaltered in 
the interaction, whilst jy gets a contribution from Sz- 



(20) 



I 

C. Canonical Operators 

In order to further emphasize the similarity in the atom 
and light evolution we will consider them on a more equal 
footing by introducing canonical variables. From the defi- 
nition of the Stokes operators for light (8)-(9) one can see 
that for a large classical field polarized in the a;-direction 
and sufficiently small angles of Faraday rotations the Sy 
and Sz components are proportional to the X and P 
quadratures of the weak quantum field in the orthogonal 
polarization (y). 

Syiz,t)^ ^{S^XL{z,t), (22) 
Sz{z,t)^ ./{S^PL{z,t). (23) 

The quadratures introduced here obey standard commu- 
tation relation [XL{z,t), Pl{z' ,t)] = i6{z — z'). 
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Next we consider canonical spin densities. Experimen- 
tally two different configurations are typically employed. 
Either the light polarization is oriented along the same 
axis as the spins of the atoms (x) or the light is orthogo- 
nal to the atomic polarization (y) . For simplicity we will 
restrict our discussion of the general method to the case 
where the light is polarized along the polarization of the 
spins. The argument we give can however easily be gener- 
alized also to the other orientations and for completeness 
we give the results for both orientations below. Let us as- 
sume that the atomic spin is polarized along the a;-axis. 
Then the commutation relation (14) averaged over the 
random positions of the atoms read 



(24) 



where {jx) — n-{j^) with n being the average linear den- 
sity of atoms. The field like canonical variables for the 
spin subsystem, 



XA{z,t) 



.hiz,t) 



3z{z,t) 



obey the canonical commutation relation 

[XA{z,t),PA{z',t)] - i5{z - z'). 



(25) 



(26) 



Now we can write the equations of motion for light 
and atoms given by (18), (21) in terms of the canoni- 
cal operators. Upon introducing the coupling constant 

At = gai ^NaNp where Np is a number of photons in 
the driving field one obtains 



d_fXL 
dz\PL 

dtyPAj 



{z,t) 
iz,t) 



Pl 




(27) 



As opposed to the initial set of equations (17), (20) these 
simplified equations are not coupled to each other, since 
we have dropped the term H^'^'^ of the interaction Hamil- 
tonian (12) due to the second rank tensor polarizability. 
This is justified in the limit of large detuning. Their 
solution can be expressed in the form of input-output 
relations for collective variables which are introduced by 



X 



%{out) 



X 



in{oui) 



'L 



dtXLiOiL),t), 
dzXAiz,0{T)). 



(28) 
(29) 



Here T is the duration of the light pulse and L is the 
length of the atomic sample. The same definitions are 
applied to the conjugated canonical variables of light and 
atoms. Finally, the input-output relations obtained from 
the equations of motion (27) read 



JDOUt 
T>OUt 



XT 

r>in 

Xa 

iDin 

J A 



P 



P. 



(30) 
(31) 



These relations constitute the foundation for the light- 
matter interface based on the Faraday rotation, see Ref. 
[6]. 

A specific example for the application of the Faraday 
interaction is the possibility to create a memory for light. 
The quantum state of a propagating pulse of light is 
thereby mapped onto the collective spin of a cloud of 
atoms. The mapping protocol applied in Ref. [16] works 
as follows: The light pulse first interacts with the atomic 
cloud such that the evolution is described by (30) and 
(31). One then measures X£"' via homodyne detection 
and uses the result for a feedback on the atomic spin. The 
feedback should substract the measurement outcome for 
j^oMt fj-Qj]2 P™* with a gain v. One can show that the 
evolution of the atoms is then described by the relations 
[161 



^A — -^A ' 

j-,lOUt _ j-,out 

^ A — ^A 



npr, 



lyXl""' = PX\l - Kiy) - lyXl 



Assuming k = = 1 wc have 



vout 
^A 



^A ' 



P, 



L I 



A — —A 



L I 



(32) 
(33) 

(34) 
(35) 



meaning that we have stored the light quadratures in 
atoms 



out^ 



{pT) 



{PD 
-{XI 



(36) 
(37) 



The quality of the mapping can be characterized by the 
fidelity T. Assuming a random set of coherent states as 
the input state the fidelity is found to be [38] 



AX' 



(38) 



Assuming light and atoms to be initially shot noise lim- 
ited we get a fidelity of J" = « 82%. The fidelity 
is limited by the initial noise of the AT™ spin component 
which was not canceled during the mapping process. If 
this quadrature would be squeezed before the pass of the 
light pulse the fidelity can approach unity. 

It is known from earlier works [6, 27] that the coupling 
constant is related to the optical depth by = drj, 
where the optical depth is d — NaCr/A, the scattering 
cross section is a = 3A^/27r and the transversal area of 
the atomic sample is A. The coefficient rj ~ Np/A'^ was 
called atomic depumping as it is closely related to a spin 
decay rate. The fact that the coupling constant depends 
on the decoherence rate is a very important observation. 
It means that for the great variety of quantum protocols 
based on such light-matter interaction a good perfor- 
mance is a compromise between the interaction strength 
and the decoherence. E.g. in order to reach a value of 
K = 1 in the memory protocol for a given optical depth 
d a nonzero atomic depumping r/ ^ 1/d is required. By 



7 



working with an ensemble with sufficiently large opti- 
cal depth d ^ 1 the atomic depumping can be made 
small and this is the basis of the protocols for quan- 
tum interfaces based on the Faraday interaction which 
have been implemented in practice. In reality however, 
there is a limit to how large an optical depth that can 
be obtained in practice. Hence there will always be some 
atomic depumping for any implementation of a quantum 
information protocol. A proper optimization of such pro- 
tocols requires knowledge of how exactly the depumping 
r] is related to the spin decay rates. In the next section 
we give the answer to this question. 



IV. DECOHERENCE 
A. Example of decoherence 

Before presenting the rigorous treatment of the general 
case of spontaneous emission for multilevel atoms, let us 
include decay in a phenomcnological basis and examine 
how an arbitrary decay would affect the atom light inter- 
action. We expect that the spins and light would decay 
according to the Heisenberg-Langevin equations 



f^L,^, ^l = {x,y}, (39) 
j^ + F,, i = {x,y,z}. (40) 



Since the decay process for light is due to scattering of 
the photons out of the mode of interest we shall sec below 
that the added Langevin noise is simpler for light than 
the one for atoms. For light the noise is just the minimal 
noise required to preserve the commutation relation for 
the field. The noise is delta correlated in time and space 
and has a vanishing mean value 



a) 



d 






2L 


d . 









(41) 
(42) 



The atomic spin decays in a number of ways. First, it can 
decay to another hyperfine level {F — 1) and disappear 
from the interaction as it is show in Fig. 4a. This is the 
same process which happens to photons. Another way for 
an atom to decay is to decay to the same hyperfine level 
as shown in Fig. 4b. This results in essentially a random 
rotation of the original spin and it creates extra noise. We 
will derive the correlators of the spin noise operators in 
the following section, as well as the corresponding decay 
rates. 

Using the definition of canonical variables for atoms 
(25) we obtain the following expressions for the decay in 



F' 

Es 

F 
F 



b) 



E, 



F' 

Es 

F 
F 



FIG. 4. Scheme of two possible decay processes. 



canonical variables 



dt 



X, 



dt 



Pa 




Jy 



1 d 



XA+Fxiz,t), 
PA + Fp{z,t). 



(43) 
(44) 



Here Tx(p) = ry(^) - ^r^,, and Fx(p) = Fy^^f / . 
Now we can consider equations of motion for the light- 
matter interaction in the presence of the decoherence in 
terms of the canonical variables 



dz\PL 



(z,t) 




iz,t), (45) 



+ (^^)(z,t). (46) 



Here the coupling constant k{z, t) is a function of position 
and time due to the decay of the classical spin component 
and of the driving wave amplitude according to (39) , (40) 



K{z,t) = ^gaiy/{j^{z,t)){a^{z,t))\^ 
- K.p-5(7=.^/i+rxt/T) 



(47) 



Integration of the equations of motion over space and 
time gives us the input-output relations 



TDOUt 

vout 

TDOUt 



XT 

Tjin 

Xf e-'^^ 



KA 



P 



A 



pin 




Fx 

Fx. 
Fp. 



(48) 
(49) 



Due to the space and time dependence of the function 
K{z,t) the integration of the equations of motion can be 
done by expanding light and atomic variables over Leg- 
endre polynomials in space and time domain [34, 39]. 
However, we consider the decay rates to be small and 
therefore take into account only the zeroth order Legen- 
dre modes for atoms and light which are just integrals 
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over length of the ensemble and over time of the interac- 
tion. The obtained relations are a generalization of (30) 
and (31) for an arbitrary decay. The following notations 
were used 



ka = K-h 



(50) 
(51) 



h{x) = -(1 -e-^). 

X 



The protocol of quantum memory based on the direct 
mapping requires k^i = 1, as described above. It is dif- 
ferent from the constraint n = \ without decoherence. 
For a given detuning an increase of the number of atoms 
leads to higher scattering for the signal field and an in- 
crease of the driving field strength gives rise to decay of 
atomic spins. Overall, it means that if the number of 
atoms for a given detuning is too low there is no way 
to meet the constraint by increasing the field power and 
vice versa. 

The noise terms in the input output relations (48) and 
(49) are a bit lengthy to express exactly. One can show 
that the noise terms have the following approximate val- 
ues 



couple to their own reservoir and do not interact with 
each other. For a justification of this approach we refer 
to Ref. [30], where it is shown that the approximation is 
suitable for dilute elongated samples, for which <C 1 
and the Fresnel number of the ensemble is much larger 
than unity. The interaction Hamiltonian including both, 
forward and non-forward modes E = Ep -|- Eg is 



Hj = E(^)ajE(+) 
— ^ttjEp'*'^ 



4 ^^Ms 



(55) 



where j labels the j-th atom. We have neglected the 



much weaker contribution E 



(-) 



(+) 



which has no en- 



hancement the by strong field in Ep. The coupling of 
forward modes to forward modes we identify as the co- 
herent interaction 



coh 



E^ 'aE 



(+) = 



A ^ 



(56) 



Here we have introduced dimensionless a which will be 
more convenient to use in what follows. The interac- 
tion with the environment Eg we will treat in Wigner- 
Weisskopf approximation. To do this we write the per- 
turbation Vj in the form 



XaI 



{F. 



XI 



6 



7y, 



{Fk)^{Fl). 



(52) 
(53) 



Where the collective noise correlators for atoms are de- 
fined as follows 



{F' 



X(P) 



1 



dzdz' I dtdt' {Fx(p){z,t)Fx(p){z' ,t')) . 



In order to preserve mean values of the canonical vari- 
ables during the mapping of the field state to the atoms 
one has to choose the feedback gain parameter to be 
V = . Assuming input states of light and atoms to be 
coherent we find the fidelity according to (38) which for 
low decay rates reads 



^(l-§7.-^[(^^l)+2(F|)-rx]). (54) 



In order to evaluate the Fidelity we need to know all 
these decay rates and Langevin noises for a particular 
atomic ensemble. Here we come to the central part of the 
paper — the proper treatment of spontaneous emission 
in multilevel atoms. 




(&L«jEp 



(+) 



E 



jbka)- 

(57) 

The summation here is performed over all directions and 
polarizations of the non-forward electromagnetic modes. 
The resulting interaction will depend on the relative ori- 
entation of the atomic dipole moments hidden in the po- 
larizability tensor a and the light field polarization vec- 
tor. The net effect of this directional dependence appears 
as a difference in the decay rates and added noises for par- 
allel and orthogonal light and atomic spin orientations. 
The modes in Eg have a Hamiltonian 



d^k 

(2^ 



uJkbLb. 



'ka- 



(58) 



We assume the system to start out in vacuum and pro- 
ceed in a Wigner-Weisskopf approach by first finding the 
time evolution of the 6-operators 



d_ 
Ft 



bkait) = i[HR + Vj,bka 
= -~iuJkbka{t) - 




^e,.a,(t)E(+)(t), 
zeo 



(59) 



which has the formal solution 



B. General method 

In this section wc find the full equations of motion 
from the single atom theory and assume that the formal- 
ism can be extended for a large collection of atoms that 



bka(t) = bkaiO)e 
l~ uJk 




dt'aj{t')-E'p^\t')e 



-iujk{t — t') 



(60) 
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We are interested in the equation of motion of some op- 
erator A (which can belong to either atoms or light). We 
insert the found expression for Es into Vj 



j^Ait) = [Hl^,A]+^'^f£ J ^^^^ka (6L(0)e^"''*KE(+\ A](t) + [Eir^a,, A](t)6,.(0)e— 

- itT^I (0g4.^'rft'(ErHt')a.(t')e^'^^-(*-*')[a.EW,A](t) - [Eir'a,, A](t)a,(t')EW(t')e-^"^^'~*')) • 

(61) 



Going into the rotating frame of the forward field, which 
has the carrier frequency ojq and performing the Markov 
approximation we may take the polarizability outside the 
integral and find that the general equation of motion for 
an observable A is given by the quantum Langevin equa- 
tion 

+ {K~'''i,Mll + tj l-A*' > -^l) ■ (62) 

Here 7 is the decay rate 7 = dp'^o/STreo which is related 
to radiative decay rate of an atom with exited state elec- 
tronic angular momentum J' as ^rad = 7/(2 J' -I- !)■ We 
have neglected the Lambshift, which is assumed to be in- 
corporated in the transition frequency ujq, also we ignore 
collective effects in the limit of a dilute ensemble (see [30] 
for details) as discussed in the beginning of Section IV B. 
The decay is described by the Lindblad form 

CM) = I (2E(-)a,Aa,E(+) 

- E^-^a^E'+^A - AE^-^af E^+') , (63) 

The decay gives two different contributions, one where 
the atom decays to the same state and will involve the a- 
coefficients belonging to app and one to different states 
through app±i and be described by ^-coefficients. But 
we are not going to have any cross terms of the form ab 
since they will oscillate at much higher frequencies cor- 
responding to the hyperfine splitting ^ GHz and average 
out to zero in (61). Above we have introduced the noise 
operators, 

(64) 

which in the Markov approximation are delta-correlated 
in time 

{[kAi),fUn)^5,iS^J{t-t'). (65) 



I 

Below we shall apply these equations of motion to ob- 
scrvables for both light and atoms. There is an essential 
difference between the decoherence of light and of atoms. 
The reason is that each time the atoms undergo sponta- 
neous emission they remain in the ensemble. For light 
on the other hand a lost photon simply attenuates the 
beam, but there is no change in the coherence of it. It 
means that we can model the resulting noise for light as 
vacuum operators when considering the field a, . 



C. Light 



In the last section we showed how to treat the spon- 
taneous emission and how to calculate the decay for the 
operators of interest. We will now apply this to the light 
field operators ax{z^ t) and ay{z^ t). First, one has to gen- 
eralize the dynamics of a single atom studied above to a 
spatially extended ensemble. The continuous equations 
of motion are obtained from the single atom version (62) 
by replacing operators with the space dependent ones de- 
fined in (13) and then integrating the right hand side over 
the ensemble volume. 

In order to obtain the required equations for the field 
propagation one has to consider the field evolution in- 
cluding the interaction Hamiltonian and the free field 
Hamiltonian. The resulting equation is obtained from the 
Heisenberg-Langevin equation (62) by replacing ^ — )■ 
■§1 + §1 ss showed in Appendix C and [40]. Since the 
light pulse is usually much longer than the size of the 
atomic ensemble we can omit the time derivative as it 
only describes retardation effect. Using the commutation 
relation for light field operators one obtains the following 
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equations of motion 
d_ 

dz JL 

_\9\lNa 

4A L 



{{a'^)^f_,af_,{z, t) + {a'^)^^a^{z, t)) 



(66) 



FL,^{z,t)^i^ ^-^^{af,f,ff,{z,t)+af,t,ft,{z,t)), (67) 

where /i — {x, y} and fi v. Here we have assumed that 
the atoms are evenly distributed so that (q;^)^^(z) = 
'7f^(«^)MM element (a^)ij is to be understood as 

the ij'th element of the matrix a^. Using (65) we obtain 
the following properties of the averaged commutators of 
the noise operators Fl.^ 



{[FLA^,t),Fljz',t')]) 



\9\lNa 2 

— ; \C( 

2A L 



{a')^,J{t-t')6{z-z'). (68) 



Since the matrix depends on the atomic operators, the 
expressions derived here incorporates that the decoher- 
ence of the light field depends on the state of the atoms. 
To evaluate the formulas we therefore need to specify the 
atomic state. Here we assume that the atoms are initially 
polarized along the x-axis and evaluate the expressions 
in a coherent spin state. In principle there is also a pos- 
sibility for cross-scattering between the x and y polar- 
ization, which is contained in the matrix element {a'^)xy 
In Appendix D we show, however, that this matrix ele- 
ment vanishes for this particular spin configuration. As 
noted above experiments sometimes use a different con- 
figuration with the atoms polarized perpendicular to the 
light polarization. Also in this case the cross polarization 
scattering vanishes, but this is not the case for a general 
orientation. Experimentally it is advantageous to avoid 
the cross polarization scattering since this avoids a rota- 
tion of the mean polarization of the light. Such rotation 
may be hard to control precisely and this classical noise 
could therefore overwhelm the quantum noise around the 
mean which we are trying to control. Therefore, the ge- 
ometry with spin polarized parallel or orthogonal to the 
direction of light polarization is preferable. In both these 
cases we obtain the equation of motion for light in the 
form (39) with the decay rates. 



7m 



|g|7 
2A 



(69) 



Faf 



Here the difference in the orientation of the atoms is again 
contained in the matrices {a^) For Cesium with F — A 
the coefficients are shown in Fig. 5. In the limit of 
large detuning the coefficients read 

A^iy) ^ 24. (70) 



25 r 




Ql -A,MHz 

500 1000 1500 2000 2500 3000 



FIG. 5. Light attenuation and noise are described by the ma- 
trix . Solid line shows the matrix element Ax as a function 
of detuning —A in MHz and the dashed line corresponds to 
Ay. We assume atoms to be polarized along x. Cs in _F = 4. 



Taking into account coherent part of the interaction 
Hamiltonian we arrive at the full equations for the light 
field 



d_ 
dz 



dz 



ay{z,t) 



2 



ax{z,t) 



aa;{z,t) [aij^_{z,t) + ia2{jx,jy}{z,t)] 



2ia2j^(2,0 ay{z,t) + y/%fy{z,t), (71) 



2VT 



ay{z,t) [aijz{z,t) - ia2{jx,jy}iz,t)] 



-2ia2]l{z,t) ax{z,t). (72) 



Since the x-polarized field is supposed to be strong and 
classical we neglect its quantum noise. 



D. Atoms 

In this section we consider atomic evolution in a similar 
way as it was done for light. To describe the decay of 
the single atomic spin we use expressions (62), (63) from 
above to obtain 



d_ 
dt 



Mz,t)^[^yc{Mz,t) 



(73) 

2aj,a]iz)E'p+\z). 

(74) 



We expect this decay to be proportional to the total flux 
of photons and we can as a good approximation ignore 
the position dependence of the light field, allowing us to 
write for the collective spin 



'^'^i^ J dz[a^j,+j,a^ - 2aj,a](z)E{+^ 



The most notable contribution comes from the x polar- 
ized part of light. Since we already expect the effect of 
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spontaneous emission to be small in the regime where the 
quantum interface is operating (remember that it scales 
as 1/A^) we can safely only consider the xx element of 
the matrix 



-jia -2ajia. 



With this assumption we find that we can make the lin- 
earization 



(75) 



where the field amplitude is given by = uJo/2ef)A 
and the decay rate coefficient is defined through the 
expression 



Sijj(i) 



dz (,{z, t). 



(76) 



This element gives the magnitude of the decay of the 
respective component of the spin. For the decay to a dif- 
ferent hyperfine level F (the h terms) , we will neglect the 
term 2ajia, which describe the increase in the popula- 
tion of the final state. The reason is that once an atom 
decays to another hyperfine level F it will no longer be 
interesting for us, since we restrict our analysis to collec- 
tive behavior of many atoms in the same F state. This 
is reasonable since the energy spacing between the two 
ground F = / ± ^ hyperfine levels is big. For the same 
reasons we obviously need to keep 2ajia for the a-terms. 
Written in short form the decay reads 



d_ 
di 



4A T 



Eiji{t) + Noise. 



(77) 



This expression gives us the decay rates of the spin com- 
ponents defined in (40). The noise terms will be treated 
in detail below. The decay rates depend on the relative 
orientation of spin and light polarizations as it is the case 
with the light decay. The expression for the decay rates 
now read 



(78) 



= 



Faf 



Similar to the discussion of spontaneous emission in sec- 
tion III C we have here expressed the atomic scattering in 
terms of the coupling constant k and the optical depth d. 
The details of the interaction and the difference between 
different atoms, which is our main interest here, is thus 
contained in the coefficient B. The required coefficients 
S are calculated in Appendix D 3 for an arbitrary atom. 
The Bi coefficients for Cesium atoms are shown in Fig. 6, 
explicit expressions can be found in Appendix E. In the 
limit of high detuning the coefficients for Cesium atoms 
with F = 4 read 

^ R n ^'^ 

37 „ „ 1 



Bx±, Bz 



:Br 




-A, MHz 



3000 



FIG. 6. Spin decay coefficients Bi as functions of detuning 
— A in MHz for parallel configuration (spin and light are x 
polarized) are shown by thick solid line for x spin component 
and by thin solid curve for y and z spin projections. The 
orthogonal configuration (spin is x polarized and the light is 
y polarized) is presented by the dashed curves, thick hne is x 
and z, thin is y spin component. The plot is calculated for 
Cesium atoms with ground state F = 4. 



The remaining term in the Heisenberg-Langevin equation 
(62) gives us the spin noise components 



Ft 



(+) 

F 



From here the noise operators for canonical spin density 
variables defined above as Fx{z, t) and Fp{z, t) are given 

by 



(79) 

These are the noise operators introduced in (46). The 
collective noise operators are obtained by averaging it 
over the light pulse duration and the length of the atomic 
ensemble 

^^(^) = ^ / / dtFxiP){z,t). (80) 

When forming the combinations such as (FxFx) then, 
since we have assumed that there are no photons in the 
reservoir with frequency at wq, the only combination 
which survives is (ff^). Moreover, we have made the 
assumption that the x (or y in case of orthogonal con- 
figuration of light and spin polarizations) component of 
the light field is dominating, so we will only have to con- 
sider the xx{yy) element of the matrix C,f = {i[a,ji\)'^. 
So ultimately what we will have left is 



{F 



X(P)/||(J_) 



r (Cf)iia) 



2 A 



p ^ "^2/(^)iia)' 



(81) 



The different (^^) are calculated in Appendix D2 and 
are generally quite complicated. Since any orientation 
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FIG. 7. Spin noise coefficients d as functions of detuning 
—A in MHz for Cesium with ground state _F = 4. The sohd 
line represents y and z noise coefficients for x polarized atoms 
and light. The dashed curves correspond to the orthogonal 
configuration of atomic spin and light. The thick dashed curve 
shows Cy^ and the thin dashed line is C^^ ■ 



besides the parallel and orthogonal will give too much 
noise on the light as discussed above, we only consider 
these two settings. The spin noise coefficients Ci required 
for the calculation of canonical spin variables noises are 
shown in Fig. 7, the explicit expressions are given in Ap- 
pendix E. In the limit of high detuning the coefficients 
for Cesium atoms with F = A read 



Ci/ii , Cz 



29 



a 



53 

y± -J^ 



37 



The cross correlations of the Langevin noise are closely 
related to the spin decay rates as shown in Appendix D 4. 
The commutator for noise reads 



([Fx,Fp])||(j_) = -(r, + r,-r,)||(j_). 



(82) 



Moreover, one can show that the anticommutator of the 
noise terms is equal to zero. 



V. SUMMARY 

In this section we summarize the obtained results and 
as an example we apply the results to determine the de- 
cay rates and Langevin noises for the specific case of Ce- 
sium atoms. First of all, the decay of the spins and light 
is described by the Heisenberg-Langevin equations 



T 



M = {x,y}, 

i = {x,y,z}. 



As it was mentioned earlier there are two preferable con- 
figurations available. Either the light polarization is ori- 
ented along the same axis as the spins of the atoms (x) 
or the light is orthogonal to the atomic polarization (y). 
In these cases the decay rates for light are given by (69) 



and the atomic decay rates can be found in (78). They 
read 



7m 



2^^A 



*IK-L) 



Fal 



From these expression we see that one can obtain a strong 
interaction k ^ 1 with negligible decoherence provided a 
that the atomic ensemble has a sufficiently high optical 
depth d ^ 1. Furthermore the decoherence of the light 
modes can be suppressed relative to the atomic decay 
provided a larger number of photons than atoms is used 
Np ^ Nji. Most of these expressions are well known 
from previous work in the field [6] . The main new result 
in this work are the exact expressions for the coefficients 
A, B and C, which are numbers of order unity, which 
contain the information specific to the particular atom 
one is considering. We will evaluate the required coeffi- 
cients for ^'^•^Cs optically pumped to ground state F ~ 4. 
In this case we have / = 7/2, J = 1/2, J' = 3/2, F = 4, 
F — 3. The coefficient oi is given by (B9), which is 
evaluated for F = F = 4 and k = 1 



ai 



7 



5760 



A45 

A 



A35 
A 



(83) 




The matrix (a^) whose diagonal components enter in 
can be found in Appendix D. The Cartesian com- 
ponents of the matrix for the spin oriented along a;-axis 
are expressed via its spherical components by (D4) and 
(D5). The required spherical components of the matrix 
are given by (D3). One needs to know the coefficients 

a^^ mentioned above and the coefficients Ck given by 
(Bll). The result for Cs atoms is 



ivv 



Since the spin is considered to be oriented along x, the 
xx-element of the scattering matrix represents decay of 
the field polarized along the spin, and the yy-component 
is correspondingly the decay rate coefficient for the or- 
thogonally polarized field configuration. 

Calculation of the spin decay coefficients S entering 
the B coefficients is similar but a bit more lengthy since 
it involves three spin projections and two possible config- 
urations. The basic expressions for parallel and orthogo- 
nal configurations of spin and light are given by (D25)- 
(D30). The required spherical components of {£,^ijv) are 
defined in (D22)-(D24). The definitions involve the same 
coefficients a^^ and cj, used above for the evaluation of 
the light decay rates. For reference we provide the result 
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of the calculation for the expression (^j,) 
atoms F = A. 



for Cesium 
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The Langevin noise operators for the light modes (68) 
are simply the minimal noise required to preserve the 
commutation relation for the free field. 

{FL,,MFl^z',t')) = ^s,jit-mz-z'). 

In the context of quantum information protocols it is 
more convenient to study the evolution of the canonical 
variables for the spin. In this case the decay process is 
described according to (43) by 

-PA = -^PA + Fp{z,t). 



dt 



T 



Here ^xiP) = ^y(z) - I'^x, and Fx{p) = Fy(^)/./{Q. 
The Langevin noise correlator for the collective canonical 
variables is found in (80) and (81) and it reads 



P d '-^W(^)|l(-L)' 



'x(p)/||(j_) - p ^ ^y(.^)i\(±)' - pg^2 

The averaged commutator for the noise is given by (82): 



{[Fx,Fp] 



IIW 



+ rp)| 



To evaluate this expression one needs to calculate the ma- 
trix {(f) . The calculation of this is essentially the same 
as for the spin decay coefficient. The expressions for the 
Cartesian components via the spherical one are given by 
(D13)-(D18). Spherical components for different F and 
F are found in (D10)-(D12). For reference we provide 
the result of the calculation for the expression {Cy)\\ for 
Cesium atoms with F = A. 
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Expressions for all of the coefficients A, B, C for ^'^■^Cs 
are explicitly given in Appendix E. 



A. Example of application 

Now we can use the derived expressions for the light 
and and atomic decay rates to evaluate the Fidelity ob- 
tained in section IV A. First, let us consider the cou- 
pling constant ka given by (51) for the case of ^'^^Cs. In 
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FIG. 8. Solution of the equation ka = 1 as a function of 
the optical depth d — Na(j/A. The thin solid line corre- 
sponds to the constraint without decoherence «; = 1. Thick 
solid and dashed curves represent respectively the solutions 
for the parallel and orthogonal atomic spin and light polariza- 
tion configurations. The plot is calculated for Cesium atoms 
with ground state _F = 4 and detuning A = —500 MHz. 



Fig. 8 we show the solution of the equation ka — ^ which 
gives us the required parameter regime for the quantum 
memory based on direct mapping protocol. One can see 
that within a wide range of values for the optical depth 
d = Nacr/A the requirement k.a — ^ coincides with the 
condition n = 1. However, when it reaches extreme val- 
ues the decay processes significantly affect the system 
evolution and one is required to provide k > 1 in order 
to fulfill the quantum memory condition for transfer of 
the the mean amplitude. 

For a given optical depth d one can express the fidelity 
(54) by the approximation of small decay probability. 
Consequently for the parameter regime where ka = ^ 
coincides with the constriction k = 1 we have 
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12 Fal ^' 
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3(Fai)2 



l-CL 



1 



CA 
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(84) 
(85) 
(86) 



In the limit of A — )■ —00 we can neglect the noise on the 
light the coefficients and the fidelity read 



CA„ 



12 

11 /41 
T V 12 



1 - CA 



1 



(87) 



In agreement with Ref. [6] the deviation of the fidelity 
from the ideal one is inversely proportional to the optical 
depth for large detunings. The results of the paper allow 
us to identify the coefficient of this proportionality. It 
shows an advantage of the orthogonal configuration of 
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FIG. 9. Optimal fidelity as a function of optical depth for Ce- 
sium with F = 4. Solid and dashed lines correspond respec- 
tively to the parallel and orthogonal configurations of atomic 
spin and light polarization. Thin solid line on top shows the 
limiting value •y/2/3. 



light and spin polarizations over the parallel one in the 
particular quantum memory protocol discussed here for 
illustration. 

For a finite detuning the fidelity can be optimized by 
varying the ratio of photon number to the number of 
atoms and by changing the detuning from the resonance. 
The resulting optimal fidelity as a function of optical 
depth is shown in Fig. 9. The corresponding optimal 
detuning and optimal ratio of photons to atoms as func- 
tions of the optical depth are shown in Fig. 10 and Fig. 11 
respectively. One has to take into account that these are 
very flat maxima, and hence the precise value of the de- 
tuning is not that important. This is illustrated in Fig. 12 
where fidelities for several optical depths are shown as a 
function of detuning. Furthermore one should also bear 
in mind that experimentally other considerations such 
as Doppler broadening may be important. Furthermore 
to illustrate our method we have for simplicity neglected 
the tensor part of the coherent interaction, and this may 
alter the conclusion reached here. 



VI. CONCLUSION 

In this paper we have found in accordance with earlier 
work that spontaneous emission occurs faster for par- 
allel polarizations than perpendicular — a consequence 
of selection rules and Clebsch-Gordan coefficients. We 
have also found the precise decay rates and noise correla- 
tions that take into account the full level structure of the 
atom and our procedure can easily be applied to similar 
systems. In this paper we have given a comprehensive 
discussion of quantum noise from spontaneous emission 
for light matter quantum interfaces based on the Fara- 
day effect. Taking into account the full level structure of 
the atoms we derive for the first time the full expressions 
for the decay and quantum noise arising from sponta- 
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FIG. 10. Optimal detuning as a function of optical depth 
for Cesium with F = 4. Solid and dashed lines correspond 
respectively to the parallel and orthogonal configurations of 
atomic spin and light polarization. 
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FIG. 11. Optimal ratio of photons to atoms as a function 
of optical depth for Cesium with F = 4. Solid and dashed 
lines correspond respectively to the parallel and orthogonal 
configurations of atomic spin and light polarization. 
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FIG. 12. Fidelity for different optical depths as function of 
detuning for Cesium with F = 4. Solid, dashed, and thin 
curves correspond to optical depths d — 30, 100, 1000 respec- 
tively. The calculation is done for the parallel configuration 
of atomic spin and light polarization. 
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neous emission. In agreement with previous treatments 
based on simplified models our theory shows that the de- 
cay and noise can be made to vanish for a sufficiently 
large optical depth of the ensemble. Given that experi- 
ments will always work with a finite optical depth it is, 
however, important to have a detailed understanding of 
the noise. In particular, in view of the rapid experimental 
advances a thorough understanding of fundamental noise 
sources will be crucial for further increasing the efficiency 
of quantum memories and for assessing the feasibility of 
new advances for light matter quantum interfaces. The 
tradeoff between the coherent dynamics and unavoidable 
fundamental losses, investigated in full detail here, will 
set the ultimate limits to the performance of future quan- 
tum networks. 
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Appendix A: Adiabatic elimination 

The projections defined in the section II will be used for 
our treatment of the operators and if we specifically apply 
it to the dipole-operator d we get (Note that we do not 
get any contributions from the terms PgdPg, PedP,, since 
d is a parity odd operator.) 

d={Pg+ Pe)d{Pg + P,) - P.dPg + PgdP, 

= dW-fd(-). (Al) 

Here d^+^ is the raising operator raising the atoms from 
the ground state \F,m) to the excited state |F',m'). In 
the following we will go into the rotating frame with re- 
spect to the laser frequency wq, so the energy of the atom 
is described by the Hamiltonian 



Inserting this into the Hamiltonian we get the expression 



(A2) 



where A^^-' is the detuning of the exited stated from the 
laser frequency. To ease notation we treat by both using 
a symbol 5 or e to represent respectively ground states 
\F,m) and exited states |F',m'). 



(A3) 



Taking the commutator with Hp^ + i/jnt we get the fol- 
lowing equation 



dt 



|5)(e| = -zAe|.g)(e| 



+ i1 



So in the adiabatic regime we neglect the time derivative 
and obtain 



e(+) . d^V 

9' 



(A5) 



OL 



(A6) 
(A7) 



Here d^^^dig^ ^ d^^Vdig-* is the dyadic vector product 
of the dipole operator with itself. 



Appendix B: Construction of the Hamiltonian 

The dipole moment d can be expanded in spherical 
components by 



d-^rf,e;-^(-l)«d_,e,. 



(Bl) 



The component of d for the F F' transition can thus 
be expressed as 

d'^+l= J2 {F'm'\d'-+^\Fm)\F'm'){Fm\el. (B2) 

q,m,m' 

The Wigner-Eckart theorem states that the matrix ele- 
ments can be expressed as 

{F'm'\dM\Fm)=Cf:^^{F'\\d\\F) 

= i-lf-'+^-'V2WTT (^^ J {F'\\d\\F). (B3) 

Some people use a different convention for the reduced 
matrix element which is related to the convention used 
here by iF\\d\\F') = V2FTl{F\\d\\F') . We have here 
used the following notation for Clebsch-Gordan coeffi- 
cients Cy^/j-,,„, = {jmj'm'\JM). It is convenient to de- 
fine new operators (Tp,p = d''p,p/{F'\\d\\F) . The ad- 
joint is defined as (Tpp, = {(T^,pY. Then the po- 
larizability tensor operator can be expressed as a = 
- EF'FFi{P\\d\\F'){F'\\d\\F)/AF'F)(Tjp'®cTpp. Using 
the above definition one obtains 



cTpp,^cTpp = {2F' + l)J2 E 



p,q n,m.m 



X I ^ ^ F' W^^ 1 F'\ \p^^p^\^^^* 

n p —m l V m q —m l ' " ' p « ^ ' 

This expression can now be split into its operator and 
tensor part by inserting the identity [32] 

EE(2fc + l)(-,w)(-,~J?)-^.A.. (B5) 

k=Ol=-k 

such that 
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r+,^ = (2i^' + l)(-l)^^+^'^(2fc + l) 



y^(-l)2f'+F-» ^•_]^^F'+2-p-m'- 



F' 1 1 1 1 -P" ^ 

-m' qm) \-q p I J \ —p m! —n 



\Fn){Fm\ 



(B6) 



The expression for the operator part can be further simphfied by evaluating the sum over p, q, and m' using the 
identity 



h k ji \ ( k h 32 \ ( h h ja] _ ( ji j2 ja \ jji 32 js 



Changing 3j-symbols to Clebsh-Gordan coefficients wc get 



a-p^, = {2F' + 1) ^-^l== ^(2A:- + 1) ^ ^ ^ [ E 



F+F' + l 



fc=l 



A; 



J Lp,q 



^3(2^+1) 

Finally, the overall polarizability can be written as a = — (2J' + l)|(J'||(i|| J)^-^ "^^a^^ (A) T^p., where 



k,F 



ar(A) = -(-ir(2fc+lK 



2i^ + l^(-l)^'(2F' + l) (J'F'I\ fJ'F'/l (f k F 



1 _ Sp, 



/A 



F J 1 ( 1 F J 1 ( 1 1 F' 1 



(B7) 



(B8) 



(B9) 



I 

The hyperfine splitting is Sp' = A — A^r'. Here we have \{ F ^ F 
used that 

~ ■- 

(F|M(+)||F')(F'||d(-)||F) = 

= (•_]^-)2(j'+/+F)+F-F' 





F' 






F' 






J 


J} 




J 





X J(2F + l)(2F + l)(2J' + l)|(J'||d||J)f. (BIO) 
The coefficients Cfc read 
Co = 1, Ci 



1 



C2 



/2F(F + 1) 
3 



(Bll) 



yiOF(F + 1)(2F - 1)(2F + 3) ' 



These coefficients are chosen in such a way that the 
resulting Hamiltonian contains spin component as irre- 
ducible tensors. The irreducible tensors are defined as 
follows 



+ (-l)'c|;;,^,,_,|Fm)(Fn|) (B13) 



Appendix C: Transformation of the light equation of 
motion 

Wc will work in the paraxial approximation assuming a 
flat transverse profile and write the forward electric field 
as 



EF(^,i) 



Wo 



Y.j^^eAaUt)e'^'^^'^'-^'> + h.c.), 

(CI) 



2eoA 

= \E\Y^e,{a,{z,t) + al{z,t)) 



^ Here |F| = ^'^^ liscvc defined the space depen 

%p = —Y^ Mlpp^C\l^iep(^el (B12) dent operators a,{z,t) = J afe^,(^)e''=-^ where it is 



" l=-k p,q 

For F = F we have 

Ml;^pp = J2Cl^:kki\Fn){Fm\ 



assumed that the different k are close to ko, so that the 
operators oscillate at the common frequency coq = \ko\. 
For the radiation field we have that 



dk 
2^ 



(C2) 
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Now we are to ready to form the EOM 



I = s — q, I' = —I, r = F + q — s, and hence 



d 

— a„{z,t) = t[Hint + Hi^,aaiz,t)], (C3) 



dk 



ikz 



dk 
2^ 



ujkakcr{t)e 



ikz 



(C4) 



But we also have from the exphcit z-dependence that 



dz 



a„iz,t) 



dk 
2tt 



ikaka{t)e"'\ 



(C5) 



which we recognize as —i[Hi^,aa{z,t)\. So we replace 
the time evolution from the radiation field by minus the 
derivative with respect to z and end up with 



d d 

^] a,a{z,t) = i[Hii-,t,a„{z,t)]. (C6) 
ot oz 



This procedure is directly applicable for Stokes opera- 
tors too. In the end we throw away the time derivative, 
which is the same as ignoring retardation effects, but this 
approimation can even be made exact by introducing a 
suitable rescaled time. By doing this we have thus an 
equation in position and not time for light observables. 



Appendix D: Calculations of atomic spin noise and 
decay rate coefficients for light and atoms 

1. Light matrix 



as 



= ^ar(A)T 



k 

FF- 



(Dl) 



k,F 



We assume that the atomic spin is polarized along x 
axis and during interaction the collective spin experiences 
merely small rotations. Having the quantization axis to 
be the x-axis we have to linearize the spin dependent 
operators around the initial state \FF). Using the ex- 
pression for the irreduccible tensors (B12) and averaging 
over the initial spin state we obtain 



FF FF 

2\ _ "fc Qfc' ^FF pFr 

2-^ r, n, , 2-^ 2-^ ^Frkl^FFk'l' 



k,k'F 



CkCk' 



2. 

p,g r,s,ld' 



(D2) 



We have changed the notation for Clebsch-Gordan co- 
efficients to the shorter form {jmj'm'\JM) = C'-^j,^,. 
Taking into account restrictions on momentum projec- 
tions given by Clebsch-Gordans the sum over four indexes 
is reduced to a sum over a single index. We have p = q, 



E 



FF FF 



FF ^F,F+q-s 

f, f,,, ^ / F ,F+q-s:k,s-q FF:k',q-s 

k,k',F " " q,s 



^ '^lq-k,s-q'^ls-k', 



g_se,®e*. (D3) 



The required elements of the matrix for the atomic 
spin aligned in the x-direction in Cartesian coordinates 
read 



= ^ + ("')++) 



Ixy 



= 0. 



(D4) 
(D5) 
(D6) 



2. Atomic Langevin noise [i[a,ji 



Let us consider the matrix C,i = i[a,ji\. The commu- 
tator of an irreducible tensor and spherical spin compo- 
nents reads 

WIff,],] = -Vk{k + ^)Clt^;: Mt+,,PF- (D7) 

This is true for an irreducible tensor defined for a single 
spin state therefore we consider the case of = i^. 



^FF 



FF 

^E7-Vfc(fcTT)E 



X Ctf+^ Ml" 



''kl^lfi 



l+li;FF'-^lp-kl'^'^ 



(D8) 



Using the same method as in the previous subsection one 
obtains 



The dimensionless tensor polarizability can be written {C^iCv 



FF 



FF„FF 



CfcCfc 



k.k' 

^k,l+fiy^k',l'+^^FF 



!^^k{k + l)k'{k' + 



E 

q r^s,l^l' 



^ ^kl;lfl ^k'l';liy'~'Fr;kd+li^FF;k'l' + 
X (^lp;kl^ls;k'l'^p'^ ^*q 



(D9) 



After some transformations one comes to the following 
expression 



{CtJ.Cv)pq 



FF FF 



J2^^^Vkik + l)k'{k' + l)J2 



k.k' 



CkCk 



c ^k,s — q—u^k',q—s+v 

X 'Jp,q+fi+v'-^k,s-p-lii'^k',q-s;lv 

y ^FF ^F,F+q-s+ii 
^ '^F,F+q-s+v;k,s-q-v'^ FF;k',q-s+v 



^ '^l,p;k,s~p'^ls;k', 



q-s 



(DIO) 



The case of F 7^ F has to be considered specifically. The 
commutator with spin components doesn't have a nice 
form anymore 

[Mt,FF,J^.] = VFiF + l) [ct,kiC^:n;i^\Fn){Fm 
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Using this expression for the commutator and doing some 
algebra we obtain 



FF FF 

^ riFF ^F,F-v \ ^ Qfc O'k' 

^ ^F,F-t,;lti^FF;lu 

k,k' 

„ riF.F+q-s+v ^F,F+q-s+v 



F,F-n-k,p-s^F,F+y;k',q-s 

p^ls:,k',q-s- (Dll) 



The full matrix (C^) is just a sum of all of these contri- 
butions for different F. 



pq 



pq 



(D12) 



Since the spin is prepared in a polarized state along the 
X direction the spin noise coefficient for i-th spin com- 
ponent in presence of the driving field polarized parallel 
to the spin is given by the xx components of the matrix 
{(f) . The case of the orthogonal configuration of the spin 
and light is described by the yy element of the matrix. 



(CoCo) 



^{(C-C+>00 + (C+C->00}: 



{(CoCo 



± — 



(CoCo)++} 

+ (C+C+>- 

{(+(-)- 

(a± = 7{(C-C->-+ + (C+C+>+- 
(C+C-> 



+ {(-(+)++ 

1 



(D13) 
(D14) 
(D15) 
(D16) 

+ (C-C+)- 

+ (C+C->++}, (D17) 

--(C-C+)- 
-(C+C-)++}. 



(D18) 



3. Spin decay rates and 



In this subsection we derive expressions for the spin 
decay rates. The coefficients for the spin components 
decays are defined via linearization of the operator — 
a^ji + jiO^ — 2ajia = — i[a, Ci] for the given polarized 
spin state. We want to have ^ '^iji which implies 

(e.)=S,0-), = S,af ). (D19) 

The last equality is always nontrivial so we can use it 
for defining Ei = Using expressions (Dl) and 

(D8) we proceed 



^FF 



-,FF„FF 



fc.fc' 



CfcCfe 



p^q r,s,l^l' 



^ \'^Fr-k'l'^Fm;k,l+ti^lp-k'l'^ls; 

CFn r^Fr y^ls /-^If? 

Fr-k,l+p,'^Fm-k'l''^lp-kl^l8;k'l' 

X \Fn){Fm\ e^® e* 



(D20) 



After averaging over the initial spin polarized state \FF) 
we obtain the required matrix elements in the spherical 
basis 



{^tijy)pq 



k,k' 



FF„FF 



k' 



-y/k{k + l)F{F + l) 



E '^ki: 



CkCk' 

l+p.f-,F,F+u 
l^l '-'FF-li' 



„ I y~tFF /^Fr ^Is f^lq 

X I Fr-k'l'^ F.F+f-k.l+p'^lp-k'l'^ls-kl 



r~iFF f^Fr f^ls i^lq 

~ '^Fr;k,l+p.^F,F+v-k'l'^lp-kl^ls;k'l' 



(D21) 



The Clebsch-Gordans require the following constraints 
on the indexes to be fulfilled for the first term r — F + 
p ~ s, I = q — s, I' — s — p and for the second term 
r — F + q — s + u, l~s — p, l' = q — s plus the usual 
condition p = q + + v. Finally, we obtain 



pq 



^p,q+p.+v 



E 

k,k' 



FF FF 



CfcCfc' 



'FF:lh 



^k{k + l)F{F + l)Y, 



C 



FF 



c 



F,F+p-s 



F,F+p-s:k' ,s-p'^ F,F+u:k,q-s+fi 



f^k,q-s+n^ls ^ig 
^ '~^k.,q-satj.^lp;k' ,s-p'-^ls-k,q- 



c 



FF 



-iF.F+q-s+l' 



F,F+q-s+y:k,s-p+^i'^ F,F+y-k\q-s 



f-ik,s~p+ti^ls 



f^lq 
0^ Is-.k' ,q—s 



(D22) 



In case of F 7^ F essentially the same calculations provide 
us with the necessary expression 



^FF 
pq 



FF FF 

sr^ at a 



"p.q+p.+f / , 
k,k' 



CkCk 



'-F{F+i)cfi+:: 



I ^) '-'lp:k,s-p'-'ls;k' 



^ \ ^F,F-s+p^F,F-s+p 



FF-k,p-s^ F,F+fj.+u-k' ,q-s^ F,F+iy;lp. 



c 



F.F+fi+u 



I f-tFF f-iF ,F+q-'S+iJ ^F.F+q-s+u 

+ ^ F,F- ti-Afi^ F,F-fi;k,p-s'~' F,F+:^;k' ,q-s 

(D23) 

The full matrix can be found by summation over 

all F 



F 



FF 
pq ■ 



(D24) 



We recall that the spin is prepared in a polarized state 
along the x direction. Then the decay rate coefficient for 
ji in presence of the driving field polarized parallel to the 
spin is given by the xx-component of the matrix {^iji). 
The case of the orthogonal configuration of the spin and 
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light is described by yy-element of the matrix. 



"^11 ~ ^(^oio)oo, 



1 

1 



{(COJO)— + (e0j0)++}: 



(D25) 

(D26) 
(D27) 
(D28) 



+ + (e+J->-- + (e+J-)++} , (D29) 

1 



2F 



{{C-J-)-+ + {C+J+)+--{C-J+) — 



{C-J+)++ - {S.+J-)— - {^+J-)++} • (D30) 



4. Relation between the spin decay and the spin 
noise 

Let us consider the commutator of the Langevin noises 
for spin polarized atoms {[Fi,Fj]). The noise is defined 
by the matrix Q = calculated in Appendix D2. 

Therefore, one can consider the mean value of the com- 
mutator {[Cy,Cz]) = -{[[<^,jy]d<^J^]])- On the over 
hand, the spin decay rate is given by the operator 

= ct'^ji + jiOt^ — '2'Oijia. In the previous subsection 
of the appendix we have shown that for spin polarized 
atoms the operator is proportional to the spin projection 
(^i) ~ 'E.i{ji). Using this fact one can show that 



{[Cy,C.]) = l^i^y + ^.-^.){j.)- (D31) 



This holds true for any cyclic permutation of x, y, z. 



Appendix E: Expressions for ^^'^Cs 



bi - 




Light coefficients 
1 



4a 



1 r 2 2 112 

2 [ao + + 56aia2 ^aoa2 



140^2 77,2 ^ /77, , 
— (6? + -62 + 2^ -6162)], 



A7, 



1 



4a 



2 L^O 



56 

ISa^ - 28aia2 + — aoa2 



539, 



9 ' ' 15 ^ V 5 
Spin decay coefficients Bf. 
1 



bib2)]. 



140 ,,2 , 77 



-ibi 



14a2ai + 49a2 
77 



6162 



1 

40!^^^! 
245 



a? - 98a2ai + 273a? 



B, 
B, 

B„ 



= B 



9 

y\\ ' 
1 



{bl + fbl+2f-^b,h 



— 2 [a^ - 14a2ai + 105a2 
4a^ 



140 



{bi 



539,2 ^ /77, , 



1 

2a 



if + 56a2ai — 35a2 



1 

175 

Is" 



/.2 693,^2 2 /77,^ ^ 
{bl + ^bl--^-b,b,)], 



(E4) 
(E5) 



4900 



(E6) 



2170 



0-2 



B^ , = B^ 



(E7) 



(E8) 

(E9) 
(ElO) 

(Ell) 



(E12) 
(E13) 



In this section we evaluate the A, B, C coefficients for 
^■^•^Cs optically pumped to the ground state F — A. In 
this case we have / = 7/2, J = 1/2, J' = 3/2, F = 4, 
F = i. In order to simplify the expressions we use the 
following notations: ak for af^ and 6fc for ■ 



ao 



ai = 



- 

144 \ 
7 

5760 



44 
21 



1 - 



A45 



3 1- 



A35 ■ 



02 



5760 



176 

16- 



A45 

A 



21 



1 - 



A45 



1 - 



A35 

A 



A35 



(El) 
(E2) 
(E3) 



Spin noise coefficients C^: 
1 



v\\ 



-2 L"i 



14a2ai + 49a2 



560 
~9~ 



77,2 „ /77, , 

—52 + 2-^— 61&2)], 



r [4a? + 308a^ 



35 

"6" 



{b\ 



1001 2 2 /77^ ,^ 



Cy 



(E14) 



(E15) 
(E16) 
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